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Abstract 

(N 

We propose a theoretical motivation to quantify actual physiological 
features, such as the shape index distributions measured by Jones and 
Palmer in cats and by Ringach in macaque monkeys. We will adopt the 
Uncertainty Principle associated to the task of detection of position and 
orientation as the main tool to provide quantitative bounds on the family 
of simple cells concretely implemented in primary visual cortex. Q 

PQ 

1 Introduction 

. 2 O ne of the fundamental tasks performed by simple cells in primary visual cor- 

,-0 tex is that of detecting position and local orientation of a stimulus [17] . On the 

other hand, the functional behavior of simple cells as visual detectors is charac- 

!— 1 terized in terms of standard linear filtering and with other so-called non classical 

behaviors [13] ■ We will concentrate on linear aspects, and consider classical re- 
ceptive profiles modeled with a planar oscillation under a spatially localizing 
H window. In [30 , such receptive profiles were studied in terms of two dimcn- 

t — > sionless indexes of shape (n x , n y ) corresponding to the product of the frequency 

of the oscillation and the sizes of the window in the direction of the oscillation 
and in the orthogonal one, showing that the distribution of such feature on VI 
simple cells of macaque monkeys is confined to a specific region. This result is 
summarized in Figure [l] Remarkably, the same confinement was found also in 
cats [2f3] |2T] , and this suggests that this pattern can be associated with some 
. . criteria of optimality with respect to perceptive tasks. Notable proposals of 

m ^ such criteria were stated in terms of sparse coding in [55] , already discussed in 

[3D] , and more recently in [5T] , or in terms of Bayesian learning [T3] . 
5^ In this paper we will focus on the task of position and orientation detec- 

tion, and propose theoretical motivations based on the Uncertainty Principle 
for the corresponding geometry to explain such confinement. In general, the 
Uncertainty Principle is indeed a tool that gives informations on the possible 
localization of functions with respect to competing symmetries, that in this case 
are those of the well known group of translations and rotations of the Euclidean 
plane. The role of symmetries in the mechanisms of visual perception in VI 
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Figure 1: Distribution of receptive profiles in terms of their shapes. Figure 
extracted from [3U] . 



is a well recognized point |5J [23] , as well as the Uncertainty Principle was 
already invoked to explain relevant cortical morphologies [10] [2] . Here we will 
use such concepts to characterize the resolution that can be obtained with joint 
spatial and angular measurements, based on the localization properties of re- 
ceptive profiles. In terms of such characterizations we will deduce the bounds 
observed in Figure [T] as the result of intrinsic notions of balance between joint 
measurements resolutions. 



2 Receptive profiles and relevant symmetries 

We will assume isotropic gaussian Gabor filters as a model for standard VI 
simple cells classical receptive profiles, defined on the Euclidean image plane: 

tfjx) = -i-e ip -^-^ e -^^ , xeR 2 (1) 

(Jy/ir 

with parameters q € K 2 , p G M 2 , a € R + . Each VI simple cell is assumed 
to perform a linear filtering with a function shaped as in ([!]), so that it can 
be characterized by these parameters. Their mapping on the two dimensional 
cortical layers are referred to as cortical maps [3J. In particular, the centers q 
of receptive fields are in a so-called retinotopic correspondence on the cortex 
[17j . while the size a is in average larger at the periphery and smaller close to 
the fovea [16] . The frequency parameters p are generally considered in polar 
coordinates p — \p\(cos6,sm8), where \p\ is called spatial frequency and the 
angle 9 up to a factor of tt is called preferred orientation, and their cortical 
maps are also well studied [25] H31 H] ■ 

The family of functions were proposed in [TUJ due to their optimal local- 
ization in space and frequency with respect to the classical Heisenberg Uncer- 
tainty Principle, and their fitness to model the linear behavior of simple cells 
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was thoroughly tested [3D]. We note however that here we are dealing with a 
simplified model of isotropic receptive fields, since as we will see this provides 
enough information for the present study, with the advantage that the results 
can be stated in a clearer form. We also recall that the real and imaginary parts 
in ([!]) correspond to so-called even and odd cells 

1 |x-g| 2 1 |jC-g| 2 

i>av( x ) = — 7= cos (p ' ( x — <i)) e 2 " 2 + i — 1= sin(p • [x — que ^ 2 

but it will be sufficient for our purposes to deal with the full complex function 
as a whole. 



2.1 Groups of transformations 

Let us introduce the following unitary operators on L 2 (M. 2 ) 

i. translations: T q f(x) = f(x — q), q £ M 2 

ii. modulations: M p f(x) = e lpx f(x), p e W 2 

iii. dilations: S ff /(x) = ), a e M+ 

iv. rotations: R a f(x) — f{r~ a x), a 6 S 2 

where r_ a stands for the usual counterclockwise rotation of an angle 9 on the 
Euclidean plane. In particular we note that rotations commute with dilations, 
and it is easy to see that 

R a M p = M TaP R a . (2) 

If we denote with 171 a L 2 (K 2 ) normalized isotropic Gaussian with unit stan- 
dard deviation 



9i(x) = 




then we can characterize the functions ([T]) in terms of the operators i., ii. and 
iii. as 

i>q,p( x ) = TgMpY.^gxix). 

Such a family is the prototype of a so-called wave packet systems [9] , and much 
is known about these structures [55] [5S] . 

In this work we will deal with the localization properties of ([!]) with respect 
to translations and local rotations, i.e. making use of the symmetries i. and iv., 
since they constitute two fundamental symmetries related to the mechanisms of 
visual perception in VI (see e.g. [8] and references therein). 

Local rotations are defined by 

Kf( x ) = f {(r q a r lx ) = T q R a T^ q f(x) (3) 

where r^x = r a (x — q) + q is a rotation of the Euclidean plane around point q, 
and with respect to these transformations we have the following. 

Lemma 2.1. Let ij)* p be as in |7p. Then 

K^l P (x) = ^ iraP (x) (4) 
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Proof. Using ^ and the definition of local rotations ^ we get 

R a^q, P { x ) = TqRaT^gTqMpE^gxix) = T q R a M p Y, a gi(x) = T q M reP R a T, a g 1 (x) 

so Q follows since rotations commute with dilations and g\ is isotropic, i.e. 
R a gi{x) = gi(x). □ 

Actually, the fact that g\ is isotropic allows to write the whole family ([I]) in 
terms of all the operators i. to iv. Indeed, denoting with 9 the polar angle ofp, 
that means p = \p\(cos9, sin. 6), we can write ([I]) as 

Tpq, P ( x ) = T q ReMnp\\E a gi(x) 

where Mi\ P \\ f(x) = e^ Xl f(x), so another way to characterize the system of 

functions |l) is to consider a family {g a> \p\ — Mn P \\'E a -gi(x) , \p\,a € R + } and 

rotate and translate each of its members. The aim of next section is actually to 
deduce properties on the localization of ip q with respect to the parameters q 
and 9, expressed in terms of the parameters \p\ and a. 



3 Measures of Uncertainty 

In this section we characterize the uncertainty associated to joint measurements 
of positions and local orientations in terms of the properties of the measurement 
devices, expressed by L 2 functions, and quantify such uncertainties for the case 
of receptive profiles. 

We recall that the generators Pj of translations along cartesian axis are given 
by partial derivatives 



1 ~ dq 



T q = d Xj , j = l,2 (5) 

q=0 



while the generator J q of a rotation around point q can be written in terms of 
the ordinary infinitesimal rotation operator J — x 2 d Xl — xid X2 



J q = — 
da 



Rl = T q JT„ q (6) 

Q— 



and acts as the skew self-adjoint operator on L (R ) 

J q f(x) = -f f {r^{x -q) + q) = ((x 2 - q 2 )d Xl - (x x - qi)3 X2 ) f{x) . 
da a=Q 

We will measure averages and variances using the standard definitions for 
operators on L 2 , denoting with (•, -)l 2 (r 2 ) the L 2 (M 2 ) scalar product and with 
|| • ||i,2(r2) the associated norm. 

Definition 3.1. Let L be a densely defined skew self-adjoint linear operator on 
T> C L 2 (M. 2 ). We define its mean value over f £T> as 

E f (L) = ((iL)fJ) L 2 (M2) el (7) 

and its variance over f € T> as 

(A f L) 2 = E f ((L - E f (L)) 2 ) = \\(L- E f (L)) /|| 2 L2(R2) . (8) 
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Since skew self-adjoint operators are the infinitesimal generators of a one 
parameter group of unitary transformations, the meaning of the average ^ is 
that of measuring the deformation of / under such transformations 

E (L) =ilim .^ expiL ^' ^ L2 ( R2 ) ~ (/</)i 2 (R 2 ) 



t-vo i 

and the imaginary constant is merely a convention to ensure the result to be real. 
With this averaging, the variance ^ has the usual meaning of strength of the 
fluctuations of / under the considered transformations, that corresponds to the 
second moment of the distribution t H> ( (exp tL) f , f) (p?) ■ This means then 
that the variance ^ provides a measure of the localization of / with respect to 
the symmetry exptL. 

When applied to the operators ^ and ^ of linear and rotational deriva- 
tives, these variances correspond respectively to a measure of linear and rota- 
tional fluctuations of a function /. The more / is insensitive to translations (/ 
smooth and close to a constant function), the smaller is its P variance, while a 
small J q variance means that / has little sensitivity to rotations around q. 

The notion of localization in orientation that arises indicates that a function 
consisting of a set of parallel stripes, independently on their widths, is maximally 
localized in orientation, while a function that is circular symmetric around q is 
minimally localized. 

If we are interested in the joint localization properties of a function with 
respect to a two parameters group of unitary transformations, generated by two 
skew self-adjoint operators L\ and L 2l we are led to consider the distribution 

(h,t 2 ) ({expt 1 L 1 )(expt 2 L 2 )f.J)L2(R2)- (9) 

In this case, if the operators L\ and L 2 do not commute, then the second 
moments of the distribution ^ are influenced by their commutator. Such an 
effect of competing symmetries is quantified by the Uncertainty Principle. 

3.1 The SE(2) Uncertainty Principle 

The operators ^ and Q satisfy the commutation relations of angular momen- 
tum [7] 

[J q ,Pi]=P 2 \ [J q ,P 2 ] = -P 2 - (10) 

These commutators define the algebra of the SE(2) group (see e.g. [5] and 
references therein) , and for them the following generalized Uncertainty Principle 



holds [HI E], with respect to the quantities of Definition 3.1 Since we are dealing 
with densely defined operators, we will skip in what follows the technicalities 
related to operator domains, and refer the statements simply to L 2 (M. 2 ). For 
more details see IT21. 



Theorem 3.2 (SE(2) Uncertainty Principle). For any f £ L 2 (R 2 ) it holds 

f {AfJ q ) (A/Pi) > \\E f {P 2 )\ 
\ (A f J«) (A f P 2 ) > WEfiP^ 



(11) 



These inequalities play the same role for the noncommutative symmetries 
of rotations and translations as the one played by the ordinary uncertainty 
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inequality for the noncommutativity of quantum mechanical operators. The 
main difference is that in this case if we consider separately each of the two 
inequalities, we can not obtain a constant lower bound. Indeed for a function / 
the product of variances of an infinitesimal rotations and a translations along one 
axis can be arbitrarily small, provided that the average of translations along the 
other axis on / is small. This effect disappears when we consider translations 
on both axis, which is natural whenever we do not want to discriminate one 
direction over the other. In this case, we can actually recast the two inequalities 



(111 into one inequality with a constant lower bound. 

The following definition is closely related to that of [I], and for this reason 
we use the same notation Angv. 

Definition 3.3. Let us define the functional 



Angv[/] 



Ef(P) 



where 



E f (P) = y/ (E f (Pi)) 2 + (E f (P 2 )y and A f P = ^{AfP.f + (A f P 2 ) 2 . 

We denote with A0[/] the corresponding measure of angular uncertainty 

AQ[f] = arctan(Angv[/]). (12) 

With this definition, a direct consequence of the SE(2) Uncertainty Principle 
is the following. 

Theorem 3.4. For all f e L 2 (M 2 ) 

(AfJ q ) Angv[/] > \. (13) 

This inequality resembles the ordinary Heisenberg uncertainty inequality, 
since the presence of a constant lower bound provides a clear constraint on the 
joint localizations quantified by AfJ q and Angv[/]. However, as first noted 



in |18j . the SE(2) uncertainty inequalities (11) can not be simultaneously min- 



imized, so also (13) does not admit minimizers. This is related to the issue 
of nonexistence of a canonically conjugate observable for angular momentum 
[TTJI12]- Indeed, if we had a well defined selfadjoint operator canonically com- 
muting with angular momentum, we would end up with a well-known complex 
equation defining minimal uncertainty states |12j , while in this case we have two 
such equations, whose solutions provide CR function functions on the R 2 x S 1 
for two noncompatible almost complex structures pQ. 



3.2 SE(2) autocorrelations 

We pass now to the study of the properties of the distribution (JoJ) applied to 
the symmetries under study, that we call autocorrelation since it has the form 
of the autocorrelation of a function with respect to the group of rotations and 
translations, and extends naturally the ordinary definition of autocorrelation 
with respect to translations. We will actually restrict the analysis to the square 
modulus of correlations, since as we will see it contains enough information for 
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the present purposes. In particular we will show that such correlations can 
be used to characterize the uncertainty in the detection of position and local 
preferred angle associated to a function. 

Definition 3.5. Given f in L 2 (R 2 ), we define its SE{2) autocorrelation cen- 
tered at q as 

C 9 [/](e,«) = |(«/,/> i2(R2) | 2 . (14) 

In general, C 9 [/](£, a) provides a natural way to study the joint localization 
properties of / with respect to position and local preferred angle. Indeed, when 
we specialize to translations we get the usual autocorrelation, and by Plancherel 
theorem 

0) = \j J(x-0W)dx\ 2 = \H\^f\ 2 m\ 2 (15) 

so we have that by Young inequality and Riemann-Lebesgue lemma C 9 [/](£, 0) is 
bounded and goes to as £ becomes large. Moreover, by the usual uncertainty 
principle we have that when / is well localized in space, then J-f is broadly 
localized, hence passing under another Fourier transform C 9 [/](£, 0) will decay 
rapidly, uniformly on q, and viceversa. 

On the other hand, if we consider correlations only with respect to rotations, 
for simplicity centered at q = 

C„[/](0,a) = I / f(r- a x)f(x)dx\ 2 (16) 



essentially the same argument applies to the decay of correlations for functions 
that are localized with respect to rotations. 

Remark 3.6 (What does "essentially the same argument" mean). 
Since TRgf — RgTf we get 



C [/](0,a) = | / f(r- a x)f(x)dx\' = | / R a {Ff){k)Ff{k)dk\' 

JR 2 JR 2 

so setting polar coordinates, with the notation </> K (<y5) = (J-f)(n cos <p, Ksintp) 
C [f](0,a)= Kdn <j> K {y-a)Mv)= ndn V e- 2mna \^(n)\ 2 

JR+ JO JR+ 



where 4> K (n) = J " e -2 ™ 1 ' 1 4> K (ip)dip and the last transition is Parseval identity. 

Since L 2 (IR 2 ) L 2 (]R + ) ® L 2 (S' 1 ) as tensor product of Hilbert spaces, and 
since f is localized with respect to rotations in the real plane if and only if it 
is localized with respect to rotations in the Fourier plane, then we can assume 
without loss of generality that </>«((/?) = f(k)^((p), where $(y) decays rapidly 
away from tp = and J R+ \r(k)\ 2 KdK = c < oo. So 

C [/](0,a) = c5>- 2 ™"* |$(n)| 2 



and now strictly the same argument used for {15) applies 
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3.3 Uncertainty associated to SE(2) measurements with 
receptive profiles 

When specialized to receptive profiles, the introduced uncertainties can be ex- 
plicitly computed. In the proofs we will use the shorthand notation g a = T* a gi, 
and 8 will be the polar angle of p. 



Lemma 3.7. The 



variance 



of the operators ^ on receptive profiles £7]) is 
Proof. Since d X] ^ a qp {x) = (ipj ~ % p {x), we get E^«(Pj) = ipj, and 



E^(P?) = \\PjW P \\h 



1 



J R 2 



y?e v dy =p 2 j 



g a (x - q)\ 2 dx 
1 

2a* ' 



□ 



Lemma 3.8. The variance of the operator |6p on receptive profiles 

\p\a 



is 



A J = Ay,. J1 = 

and we will call it angular momentum variance. 
Proof. Since 



V2 



(17) 



then 



J 9 ipqp{ x ) = i (( x 2 - qi)pi - {xi - qi)Pi) % p {x) 



Its mean value vanishes on ip° , due to the isotropy of g a : 

E 4Z v {J q ) = (^VooiVOlW) =* / (X2P1 ~ xipz) \g a {x)\ 2 dx 

Jr 2 

= i\p\ (r-ex) 2 \g a (x)\ 2 dx = i\p\ / x 2 \g a {x)\ 2 dx = 



7 .2 



M 2 



To compute the variance, by analogous arguments 



(^WT = -W ((r- e x) 2 y\gAx)\ 2 dx 



\P\ 2 

\P?Q 2 
1,Fk 



e ° 2 dx\ 

2 

y{-2ye~ v )dy 



2 — — i 

x 2 e " 2 dx2 

2 



\p\a 
V2 



□ 
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We have then obtained the following proposition, which shows that for re- 
ceptive profiles the angular momentum variance AJ is inversely proportional 
to the angular uncertainty quantified in terms of Angv. 

Proposition 3.9. Let tp^ p be as in £/p. Then 

1 

7i ' 



Proof. Using Definition |3.3[ Lemma 3/7 and Lemma 3J3 we have 

11 

<r\p\ ~ TTaj 



<Angv) Wp ] - ^ = . (18) 

□ 



We will now consider SE(2) autocorrelations of receptive profiles, and see 
that they indeed contain precisely the desired joint information on localizations 
in space and local orientation associated to the uncertainties we computed. 

Proposition 3.10. Let i\) a qv be defined by |7J). Then its SE (2) -autocorrelation 
reads 

C q W QP ]{La) = eSe-^Q--™ «>. (19) 



Proof. By Lemma 2.1 and computing the Fourier transform of a gaussian 

1 



e 



— s -«(p-(roP))x e IP^—dx 



= e 2 ? e 4° 2 e 4 

so the result follows since \p — (r a p)\ 2 — 2\p\ 2 (l — cos a). □ 

This proposition shows that the decay of the autocorrelation in space is 
a Gaussian with the same width of the corresponding receptive profile, that 
characterize spatial uncertainty. With respect to rotations, we have ended up 
with a Von Mises distribution in orientations. Such distributions appear nat- 
urally when discussing the the SE(2) uncertainty principle [7J I15j . but they 
also provide a good model for orientation tuning of simple cells |32j . which is 
defined as the response curve of a cell to oriented stimuli [32 [27] . This confirms 
that the introduced notion of localization is compatible with the resolution of 
measurements performed with receptive profiles. Moreover, we note that the 
commonly used circular variance [19; of the Von Mises distribution in (191, up 
to a normalization constant, is 

CircVar(A^) = 1 h « AJ ^ 



which results to be numerically close to what we have introduced as angular 
uncertainty ( 12 ) when applied to receptive profiles ( 18 ) 

( ae wj) 2 =H°(vikr 2 
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In particular, as we will see in next section, typical values of AJ in the fil- 
ters encountered in VI are around 1.7, where the difference between these two 
notions of variance is around 5 • 10~ 2 . 



4 Bounds on the shape index induced by Un- 
certainty 

In this section we will use the measures of uncertainty referred to receptive 



A receptive profile ip^ p consists of an oscillation of frequency v = l^j under 
a gaussian bell of width a, so it appears natural to define a dimensionless index 



profiles ( 18 1 and ( 19 1 to deduce relevant features about the physiological data 
measured in |30j and depicted in Figure. In particular we will see how the 
informations provided by the analysis of uncertainty relations of Section [3] are 
sufficient to establish bounds on the number of subregions observed in the family 
of filters implemented in VI, and permit to reobtain characteristic sampling 
rates commonly used in image analysis. 

A reci 
a gaussia 
of shape [3D] 

n = va. (20) 

This quantity is related to the number of subregions defining a receptive pro- 
files, since if we let be the number of half wavelength of receptive profile's 
oscillation within k standard deviations a, we obtain — Akn. As it is appar- 
ent from the data measured in [30j . we see that approximately k = 2 standard 
deviations are sufficient to represent the main content of the filters, so that we 
can relate the effective subregions N to n as N = 8n, 



In terms of n, the angular momentum variance (17) of a receptive profile 
ip qp reads 

AJ = V2nn (21) 
while its angular variance (12), after (18), reads 



AB[r,;,]=a,r| i m(^) . (22) 

4.1 Lower bound for orientation measurements 

By the discussions in Section [3j we have seen how we can quantify with AG the 
angle resolution allowed by a linear filtering. If we refer to the task of orientation 
detection, we can set as a reasonable bound that of angle uncertainty less than 
7r/2, that is expressed by 

Aeh// 7 1 < -. 



This condition can be stated in terms of the shape index using ( 22 ) 



1 

n > — « 0.16 = n min . 

As we can see in Figure and by the discussions in [30], cells which show a 
selectivity in orientation all lie above this threshold. Moreover, we note that for 
indexes n < n m i n it can be a hard task to distinguish an even cell from being 
represented only by a gaussian, while odd cells under this threshold all appear 
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identical up to a multiplicative factor, so the parametric fit of the Gabor model 
([!]) is quite delicate in this region. We can then interpret the bunch of broadly 
tuned cells around the zero value of the shape index n as generally below the 
minimal uncertainty bound that allows a consistent detection of orientations. 



4.2 Upper bound 

In order to discuss the upper bound, we introduce a notion of characteristic 
length associated to a specific level set of the correlations 0, intrinsically related 
to the task of detection of positions and local orientations. Its purpose is to 
quantify the minimum distance that one needs to cover in order to decorrelate a 
function / as much as / is decorrelated when compared at orthogonal directions. 

Definition 4.1. The correlation length for f G L 2 (R 2 ) is the smallest distance 
A for which 

C[/](£,0)<e[/](0,|) V|£| = A. (23) 
If we apply this notion to receptive profiles ([lj we obtain the following. 



Proposition 4.2. The shape index \20\ is bounded from above by the ratio of 
the correlation length A and the spatial uncertainty a 

2im < - . (24) 
a 



Proof. Condition (23 1 on receptive profiles tp^ p , by (19 1 reads A J < 2 ct> smcc 



so (24) follows by the relation (21) between AJ and the shape index. □ 



On the other hand, as discussed when dealing with the relation between the 
shape index and the number of subregions, we have also that the effective field 
of influence of a receptive profile can be set within two standard deviations a. 
From this point of view, we can then assume that the distance d at which a 
receptive profile is effectively spatially uncorrelated corresponds to the distance 
that one has to cover in order to let its effective effective field of influence not 
intersect with its translation at a distance d, i.e. d = Act. 

In order to couple with both position and orientation measurements, we 
will then consider the hypothesis of balance of the two characteristic scales 



introduced, that is the identification A = d. By (24), this condition can be 
stated in terms of the shape index as 

2 

n < - w 0.64 = n max . 

7T 

To compare this bound with Figure, we recall that here we are dealing with 
the simplified model of isotropic receptive fields, while in |30) the analysis is 
performed considering two anisotropic indexes n x and n y . In terms of such 
indexes we can see that the largest part of the population lies within two bounds 
n x < 0.5 and n y < 0.76, and n max looks in good accordance with their mean 
value. 
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The question of whether this identification of characteristic distances is truly 
implemented in the cortex cannot be answered at this point, but we note that 
a cortical scale related to the symmetries under study that is possibly com- 
patible with the proposed relation is the mean correlation length of orientation 
preference maps (see e.g. [21 [2] and references therein). Indeed, by the mea- 
surements performed in [5] we see that such scale is comparable with the size 
of a so-called cortical point image, that is the cortical region that is activated 
after a highly spatially localized stimulus, and at least when we reduce to linear 
behavior of cells this notion corresponds to what we have indicated as effective 
field of influence. 

4.3 Sampling on orientations 

Another intriguing consequence of the performed uncertainty analysis can be 
stated in terms of optimal sampling rates for orientation detection. Indeed if 
we consider the mean value on shape index measured in [30] , or equivalently, in 
terms of the deduced bounds, for n — " max +" mi ° sa 0.4, we have that 



With respect to Gabor filters possessing such n, one way to use such result 
is to consider that the detection of orientations at angles that are closer than 
this uncertainty do not provide an actual improvement in the resolution of the 
local orientation present in the stimulus, so that it can be sufficient to cover 
the interval of orientations [0, 7r) with a sampling having a ir/8 spacing. This 
actually compares well with the notions of optimal sampling adopted in image 
analysis tasks (see e.g. in [26] and references therein), generally justified with 
independent arguments. Moreover, this uncertainty analysis permits to set clear 
sampling spacings depending on the shape index of the filter used. 

5 Conclusions 

In this paper we have studied theoretical aspects of an analytic characterization 
of uncertainty that generalizes the well known Heisenberg Uncertainty Principle 
to the symmetries associated with the task of joint measurements of position and 
local orientation. The implications of this analysis, together with an hypothesis 
of balance between characteristic correlation distances, allowed us to obtain 
bounds comparable with experimental data on the shape index of the VI simple 
cells that are selective for orientation, and to separate them from broadly tuned 
cells, which lie below the uncertainty bound for consistent orientation detection. 

We remark that this was possible even if our working assumptions on the 
functional behavior of simple cells were reduced to linear filtering with symmet- 
ric receptive fields, and the only considered task is the one associated to the sole 
symmetries of rotations and translations. 

Whether such elementary principles could be directly responsible of the ob- 
served distribution of receptive profiles is a question that can hardly find an 
answer. Nevertheless, the present study shows that they are sufficient to de- 
scribe many of the relevant features that concern the shape of simple cells. 




12 



References 

D. Barbieri, G. Citti, Coherent states of the Euclidean Motion Group and 



CR regularity. Submitted, http://arxiv.org/abs/1301.3783 



D. Barbieri, G. Citti, G. Sanguinetti, A. Sarti, An uncertainty principle 
underlying the functional architecture of VI. J. Physiol. Paris 106(5-6):183- 
193 (2012). 

A. Basole, L. E. White, D. Fitzpatrick, Mapping multiple features in the 
population response of visual cortex. Nature 423:986-990 (2003). 

E. Breitenberger, Uncertainty measures and uncertainty relations for angle 
observables. Found. Phys. 15:353-364 (1985). 

P. Chossat, O. Faugeras, Hyperbolic planforms in relation to visual edges 
and textures perception. PLoS Comput. Biol. 5(12):1-16 (2009). 

W. H. Bosking, J. C. Crowley, D. Fitzpatrick, Spatial coding of position and 
orientation in primary visual cortex. Nature Neurosc. 5:874-882 (2009). 

P. Carruthers, M. M. Nieto, Phase and angle variables in quantum mechan- 
ics. Rev. Mod. Phys. 40:441-440 (1968). 

G. Citti, A. Sarti, A cortical based model of perceptual completion in the 
roto-translation space. J. Math. Imag. Vis. 24(3):307-326 (2006) 

A. Cordoba, C. Feffcrman, Wave packets and Fourier integral operators. 
Comm. Partial Diff. Eq. 3:979-1005 (1978). 

J. G. Daugman, Uncertainty relation for resolution in space, spatial fre- 
quency, and orientation optimized by two dimensional visual cortical filters. 
J. Opt. Soc. Am. A 2:1160-1169 (1985). 

D. A. Dubin, M. A. Hennings, T. B. Smith, Quantization in polar coordi- 
nates and the phase operator. Publ. RIMS Kyoto 30:479-532 (1994). 

G. B. Folland, A. Sitaram, The uncertainty principle: a mathematical sur- 
vey. J. Fourier Anal. Appl. 3:207-238. (1997) 

N. V. Graham, Beyond multiple pattern analyzers modeled as linear filters 
(as classical VI simple cells): Useful additions of the last 25 years. Vision 
Research 51:13971430 (2011). 

H. Hosoya, Multinomial bayesian learning for modeling classical and non- 
classical receptive field properties. Neural Comput. 24(8):2119-2150 (2012). 

Z. Hradil, J. Rehacek, Z. Bouchal, R. Celechovsky, L. L. Sanchez-Coto, 
Minimum uncertainty measurement of angle and angular momentum. Phys. 
Rev. Lett. 97:243601 (2006). 

D.H. Hubel, T.N. Wiesel, Uniformity of monkey striate cortex: a parallel 
relationship between field size, scatter, and magnification factor. J. Comp. 
Neurol. 158(3):295-305 (1974). 



13 



[17] D.H. Hubel, T.N. Wiesel, Functional architecture of macaque monkey visual 
cortex. Proc. R. Soc. Lond. B. 198:1-59 (1977). 

[18] R. Jackiw, Minimum Uncertainty Product, Number Phase Uncertainty 
Product and Coherent States. J. Math. Phys. 9:339 (1968). 

[19] S. R. Jammalamadaka, A. Scngupta, Topics in circular statistics. World 
Scientific (2001). 

[20] J. P. Jones, L. A. Palmer, An evaluation of the two-dimensional Gabor 
filter model of simple receptive fields in cat striate cortex. J Neurophysiol 
58:12331258 (1987). 

[21] J. P. Jones, L. A. Palmer, The two-dimensional spatial structure of simple 
receptive fields in cat striate cortex. J Neurophysiol 58:11871211 (1987). 

[22] C. Kalisa, B. Torresani, N-dimensional affine Weyl-Heisenberg wavelets. 
Ann. Inst. H. Poincare Phys. Theor., 59(2):201236 (1993). 

[23] H. A. Kastrup, Quantization of the canonically conjugate pair angle and 
orbital angular momentum. Phys. Rev. A 73:052104 (2006). 

[24] W. Keil, F. Wolf, Coverage, continuity and visual cortical architecture. Neu- 
ral Systems and Circuits l(17):doi:10.1186/2042-1001-l-17 (2011). 

[25] D. Labate, G. Weiss, E. Wilson, An approach to the study of wave packet 
systems. Contemp. Math., Wavelets, Frames Operator Theory 345:215235 
(2004). 

[26] T. S. Lee, Image representation using 2D Gabor Wavelets. IEEE Transac- 
tions on Pattern Analysis and Machine Intelligence 18:959-971(1996). 

[27] F. Mooser, W. H. Bosking, D. Fitzpatrick, A morphological basis for orien- 
tation tuning in primary visual cortex. Nature Neurosc. 7:872-879 (2004). 

[28] K. Ohki, S. Chung, P. Karal, M. Hubener, T. Bonhocffcr, R. C. Rcid, 
Highly ordered arrangement of single neurons in orientation pinwheels. Na- 
ture 442:925-928 (2006). 

[29] B. A. Olshausen, D. J., Field, Sparse coding with an overcomplete basis set: 
a strategy employed by VI?. Vision Res 37:33113325 (1997). 

[30] D. L. Ringach, Spatial structure and symmetry of simple cell receptive fields 
in Macaque primary visual cortex. J. Neurophysiol. 88:455-463 (2002). 

[31] J. Shelton, P. Sterne, J. Bornschcin, A. S. Sheikh, J. Luecke, Why MCA? 
Nonlinear sparse coding with spike- and- slab prior for neurally plausible im- 
age encoding. Adv. Neural Inf. Proc. Systems 25:2285-2293 (2012). 

[32] N. V. Swindalc, Orientation tuning curves: empirical description and esti- 
mation of parameters. Biol. Cybern. 78:4556 (1998) 



Davide Barbieri: CAMS, EHESS, Paris, davide.barbieri8@gmail.com 
GlOVANNA Citti: Dept. Mathematics, Bologna, giovanna.citti@unibo.it 
Alessandro Sarti: CAMS, EHESS, Paris, alessandro.sarti@ehess.fr 



14 



